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Abstract 

We consider internal diffusion limited aggregation in dimension larger than or equal 
to two. This is a random cluster growth model, where random walks start at the origin 
of the d-dimensional lattice, one at a time, and stop moving when reaching a site that 
is not occupied by previous walks. It is known that the asymptotic shape of the cluster 
is a sphere. When dimension is two or more, we have shown in a previous paper 
that the inner (resp. outer) fluctuations of its radius is at most of order log(radius) 
(resp. log 2 (radius)). Using the same approach, we improve the upper bound on the 
inner fluctuation to \J log (radius) when d is larger than or equal to three. The inner 
fluctuation is then used to obtain a similar upper bound on the outer fluctuation. 
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1 Introduction 

This note is a companion to our paper pQ. There, we introduced a family of cluster growth 
models with a spherical asymptotic shape but a wide diversity of shape fluctuations. Internal 
diffusion limited aggregation (internal DLA) was one member of this family. More precisely, 
the internal DLA cluster of volume N, say A(N), is obtained inductively as follows. Initially, 
we assume that the explored region is empty, that is -A(O) = 0. Then, consider N independent 
discrete-time random walks Si, ... , Sjv starting from 0. Assume A(k — 1) is obtained, and 
define 

r k = mi{t>0: S k (t)gA(k-l)}, and A(k) = A(k - 1) U {S k (r k )}. (1.1) 

We call explorers the random walks obeying the aggregation rule (jl.ip . We say that the k-th 
explorer is settled on S k (r k ) after time r k , and is unsettled before time r k . The cluster A(N) 
is interpreted as the positions of the N settled explorers. 

In this paper we show how the tools developed in pQ lead in dimension d > 3 to sharper 
estimates on the fluctuations of A(N) with respect to its spherical asymptotic shape. We 
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keep the notation of [T], and recall the basic ones to make the paper as self-contained as 
possible. We denote with || • || the euclidean norm on M. d . For any x in IR d and r in R, set 

B(x,r) = {y G M. d : \\y - x\\ < r} and M(x,r) = B(x,r) f]Z d . (1.2) 

For A C 7j d , |A| denotes the number of sites in A, and the boundary of A is dA = {z (jL A : 
3y G A, \\y - z\\ = 1}. For a simple random walk, let H(A) denotes its first hitting time of 
A. The inner error 5j(n) is such that 

n - 6i(n) = sup {r > : 1(0, r) c A(|B(0, n)\)} . (1.3) 
Also, the outer error 8o{n) is such that 

n + 5 (n) =inf {r > : A(|1(0, n)|) C 1(0, r)} . (1.4) 
Our main result is as follows. 

Proposition 1.1 There are constants {ad, (3d, d > 3} such that in dimension d > 3, with 
probability 1, 

hmsup < ad, ana hmsup < Pd- (1-5) 

\/log(n) ylog(n) 

Remark 1.2 For d = 2 we show, with similar computations, that there are constants a^, 02 
such that, with probability 1, 

hmsup - — -— < «2, and hmsup - — -— < p2- (1-u) 
log(ra) log(n) 

The inner error bound in (11.61) was already obtained in all dimensions in pQ. Recently, 
Jerison, Levine, and Sheffield [2] have established, in dimension two and with a different 
method, the estimates (jl.6p . Also, they announced in [2] that the approach they followed 
could be adapted in dimension d > 3 to get (II. 5p . 

Let us describe the main steps. The inner error is at the heart of the argument. It is 
based on a large deviation estimate which refines our previous estimates, with interest of its 
own. For a real x, let [^J be the integer part of x. 

Lemma 1.3 Choose R large enough, and A > 1. Assume that \_AR d \ explorers lie initially 
on 1(0, R/2). We call rj the initial configuration of these explorers and A(n) the cluster they 
produce. There are positive constants {k^, d > 2} independent of R and A, such that when 
d>3 

P (1(0, R) <jt A{r,)) < exp (-K d AR 2 ) , (1.7) 

and when d = 2, we have 

R) £ A{ V )) < exp (-Ka^) • (1-8) 
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Remark 1.4 The reason behind the previous Lemma, in d > 3, is that out of |_^4-R d J 
explorers, only about AR 2 eventually hit a fixed site on the boundary of M(0,R), so that 
it is only these very explorers that need to be pushed away from this very site. The cost 
should be proportional to AR 2 . 

For the outer error, we use a large deviation estimate symmetrical to Lemma 11.31 as 
well as our coupling between internal DLA and the flashing process of [T]. This last large 
deviation estimate was recently proved by Jerison, Levine and Sheffield in [2]. 

Lemma 1.5 [Lemma A of Jerison et al.fE/] For (3, and R positive reals, assume that \_/3R d \ 
explorers lie initially outside M(0,R). We call n the initial configuration of these explorers 
and A{rf) the cluster they produce. There are positive constants {n' d , d > 2}, such that for 
(3 small enough, we have when d> 3, 

P(0eA( V ))<exp(-K' d R 2 ), (1.9) 

whereas when d = 2, we have 

P(0E A( V )) < exp L^J^j . (1.10) 

We give an alternative proof of this result, based on estimating the probability of crossing a 
shell, while avoiding traps. 

Lemma 1.6 Consider d > 2. Fix a positive real R, and start a random walk on z e 
dM(0,2R). There are positive constants {n d , a d } such that for any V subset of the shell 
S = 1(0, 2R)\M(0, R), we have 

P z (H(B(0,R)) <H(V C )) < exp (a d - k^J where p 6 - 1 = M. (1.11) 

Remark 1.7 V c = S\V is interpreted as traps. Note that p is proportional to the radius 
of a cylinder of height R and volume \V\. We can also read (jl.lip in the following way. 

P z (#(B(0, R)) < H{V C )) < exp (^a d - n d (J^J " " j . (1.12) 

This shows that for (11.121) to be an effective inequality, one needs that \V\ be smaller than 
R d . The power l/(d — 1) on R d /\V\ in (I1.12p is not important in proving Lemma 11.51 If 
one were willing to accept the weaker power 1/d, then one would have the following simple 
heuristics in dimension d > 3. Let t denote the time the walk spends in the annulus of 
height R. On one hand, the central limit scaling yields that this probability of such a stay 
is of order exp(—cR 2 /t). On the other hand, all this time should be spent on sites of V, and 
it is well known that the probability is of order exp(— Kdt/\V\ 2 ^ d ). Putting together these 
opposite requirements, and optimizing over t, we find a statement weaker than (11.121) . but 
sufficient for our present purpose. 

P z (H(B(0, R)) < H(V C )) < exp | a d - K d (J^ ) . (1.13) 
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Even though it is not written in pQ, inequality (11.13P was the motivation behind the intro- 
duction of flashing processes in [I], which were basically used to bypass this type of estimate. 
In this paper we show how the use of flashing explorers leads easily to Lemma 11.61 

The rest of the paper is organized as follows. In Section [2], we enounce some known results: 
we recall the approach of Lawler, Bramson and Griffeath [3J, and useful large deviation 
estimates. Then, the inner error estimate is proved in Section |3j In Section HI we show how 
a flashing process permits a simple control on the outer error. Finally, we have gathered in 
an Appendix the proof of the large deviations Lemmas 11.31 11.51 and 11.61 

2 Prerequisites 
2.1 Notation 

We recall some notation of pQ. The state space of configurations is N zd , its elements are 
denoted r\ and represent starting conditions for a set of explorers, or random walks. Two 
types of initial configurations play an important role here: (i) the configuration nl z * formed 
by n trajectories starting on a given site z*, (ii) for A C Z d , the configuration 1a that we 
simply identify with A. For any configuration rj G N z<i we write 



Definition 2.1 Let R G 1R+ U {oo}. For z G M(0,R) U dE(0,R), we denote by M R (r),z) 
(resp. Wr(t), z)) the number of simple random walks (resp. explorers) initially on rj that hit 
z when or before exiting B(0, R). Thus, when z G <9B(0, R), M R (r), z) (resp. W R (r], z)) is the 
number of simple random walks (resp. explorers) which exit B(0, R) exactly on z. 

Remark 2.2 Note that trajectories of walkers and explorers can be coupled to be the same 
up to the settling time of the explorer, the walker then proceeding along its simple random 
walk trajectories. 

As in [5] (Section 3), it is useful to stop explorers as they reach dB(0,R), for some R > 0, 
then to define A R {rj) as the set of positions of settled explorers. 

Definition 2.3 Consider R G R U {oo}. We set 




(2.1) 



z£Z d 



Wz G B(0,i?), 



M R ( V , z) = W R (r], z) + M R (A R (r)), z). 



(2.2) 



Finally, for any function F : Z d — > R, and A C Z d , 



we denote 



F(A) = 
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2.2 On a classical approach 

We recall the approach of Lawler, Bramson and Griffeath in [3] . Send N — |B(0, n) | explorers 
from the origin. The approach of [3] is based on the following observations, (i) If explorers did 
not settle, they would just be independent random walks; (ii) exactly one explorer occupies 
each site of the cluster. Then, observations (i) and (ii) imply that for any integer n and 
z G 1(0, n) 

~ law 

M n (N 1 , z) := W n (N 1 , z) + M n (A n (N), z) > M n (N 1 , z). (2.3) 
When z G <9B(0, n) inequality (12.31) becomes an equality 

W n (N I , z) + M n (A n (N), z) = M n (N 1 , *). (2.4) 

Note that for any set A C B(0, n), M„(A, z) is a sum of independent Bernoulli variables. 
Note also that A„(iV) C 1(0, n) so that for any z G 1(0, n) U 358(0, n) 

W n {N 1 , + M n (l(0, n), z) > M n {N I , «). (2.5) 

However, Lawler et al. did not use that W n (Af lo,^) and M n (l(0, n), z) were independent. 
They could only obtain a rough estimate on the lower tail of W n (N lo, z). This in turn gave 
some estimates on the inner error, which was used to derive bounds on the outer error, by 
using that the cluster covers 1(0, n — 5i(n)). In other words, from (|2.4p . and the definition 
of 5j(n), for R > n and z G 31(0, R) 

W R (N lo, + Mr(B(0, n - 6j(n)), z) < M R (N 1 , z). (2.6) 

Therefore, if 5/(n) is likely to be smaller than r < n < R, and 2 G 3B(0, R), we have 

l{5r(n)<r} (W R (N I , + M fl (B(0, n-r),z)< M R {N 1 , z)) . (2.7) 

We will also make use of the independence of the a- fields generated by the events {Sj(n) < r} 
and the random variables W R (N lo, z ) on the one hand, and that generated by the random 
variables Mr(B(0, n — r),z) on the other. 



2.3 On sums of Bernoulli variables 

Let us now recall a simple tool of pQ in estimating deviations in view of (12. 5p . and (12. 6p . We 
first enounce the lower tail estimate. 

Lemma 2.4 Suppose that a sequence of random variables {W n ,M n ,L n ,M n n G N}, and a 
sequence of real numbers {c n , n G N}, satisfy for each nGN 

W n + L n + c n > M n , and M n = M n . (2.8) 
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Assume that W n and L n are independent, and that L n and M n both are sums of independent 
Bernoulli variables. Assume that the Bernoulli variables {Yj , . . . , Y^} summing up to L n , 
satisfy for some k > 1 



(HI ) sup sup E [K/ n) ] < - — - , 

n i<N n K 

(H2) fi n := E[M n ] - E[L n ] > 0. 

Then, for any n inN and £ n in M. we have for all A > 

/ N ^ 



P {W n < Cn) < exp ( -A(/i n - £„ - C„) + y f /i„ + K ^[Y^f \ \ . (2.9) 

The upper tail estimate needs other assumptions. 

Lemma 2.5 Assume for each neN, and /or an event A n 

X An (W n + L n )<M n , and M n = M n . (2.10) 

Assume that W n and L n are independent, l An and L n are independent, and that L n and 
M n both are sums of independent Bernoulli variables such that \x n := E[M n ] — E[L n ] > 0. 
Then, for all n in N, £ n in K and A G [0, log 2], 

P (W n > a, An) < exp ^-A(£ n - Mn) + A 2 ^„ + 4 ^ ^K (n) ] 2 j j . (2.11) 

Remark 2.6 This lower (resp. upper) tail estimate turns out to be useful when £ n + c n is 
less than (resp. £ n is more than) i?[M n ] — E[L n }. By Lemmas 12.41 and 12.51 tail estimates 
reduce to a three step strategy: i) estimation of E[M n ] — E[L n ] ii) estimation of E 2 \Y^] 
hi) optimization in A. We emphasize that, in particular for the lower tail, this strategy does 
not require any control of the variance of W n . 

Proof of Lemmas 12.41 and 12.5b As in [TJ this is an application of Lemma 2.3 of pQ. For 

the lower tail, using the exponential Chebyshev's inequality, then the independence between 
W n and L n , formula (12. 8p and centering the random variables, we get 

P(W n < O < ^ [e _ A( ^_ £[ ,„ 1)] J e- A(E[M " 1 - g ^^-^. (2.12) 
With, for all t G R, f(t) = e* - (1 +t) and g(t) = (e* - l) 2 , by Lemma 2.3 of P, 

E \ e -\{M n -E[M n m f JVn 

^-a^-^])] < ex P |/(-A)(^[M B ] - + ^(-A) J> 2 [lf 

We conclude by observing that for all t G 1, 



t 2 



f(t) < L. e W+ an d g (t) < t 2 e 2 ^+, (2.14) 
where [■] + stands for the positive part. The proof for the upper tail is similar. 
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2.4 On a discrete mean value property of Green's function 



Proposition 2.7 Consider d > 2. There is a constant Kd such that, for any n and R with 
n — n3 < R < n and z in 1(0, R) with n — \\z\\ < 1, 



l(0,R)\G n (0,z)- \ G n (y,z) 



(2.15) 



Proof : For n — R large enough (larger than some constant that depends only on d) this is 
Theorem 5.2 of pQ. For n = R this is a direct consequence of Lemmas 2 and 3 of [5]. For 
the remaining cases one can use the same Lemmas in conjunction with Lemma 5 of [5]. 

Remark 2.8 For the inner bound we will use Proposition 12.71 with R = n. For the outer 
bound we will use Proposition 12.71 with n — R of order log n in dimension 2 and A/log n in 
dimension d > 3. 



3 Inner error 

3.1 Exploration by waves 

We choose the following height sequence. For any positive integer n, h(n) = ^\og(n) in 
d > 3, and h(n) = log(n) in d = 2. We partition Z d into concentric shells of heights h(n). 
We define S = 1(0, h(n)), and for k > 1, 

S k = M(0,(k + l)h(n))\M(0,kh(n)), and S fc = 91(0, kh(n)). (3.1) 

We realize the internal DLA with N = |1(0, n)| explorers as an exploration wave process, 
where concentric shells are covered in turn (see Section 3 of [5]). 

We fix an integer k. For a site z e we call ce/Z centered on z, C(z) := l(z, /i(n)) n«Sfe, 
and we call tile centered on z, T(z) := l(z, h(n)/2) fl A generic cell is denoted C, and a 
generic tile is denoted T. Note the obvious facts 

|J M{z,h{n)) D<S fc . (3.2) 

Before covering shell «Sfc, one stops the unsettled explorers on E^. Following PQ, for z G Sfe, 
we prove that the W / fch(n)(-^ !o,T) explorers stopped on T = T(-z) are likely to cover C(z), 
if kh(n) < n — Ah(n) for a large enough constant A. More precisely, we show that the 
probability of the event {S k <£. A(N)} is smaller, for A large enough, than any given power 
of n. As first observed in [3] 

W kh(n) (N I , T) + M M(n) (l(0, kh(n)), T) > M kh{n) {N 1 , T). (3.3) 

Since (13.31) corresponds to an inequality of type ( 12.81) . we wish to use Lemma [231 but we 
need to ensure (HI) and (H2). 
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First, if B(r) denotes the sites of B(0, kh(n)) at a distance less than r from 7~, there is L 
and pd > 1, (which depend only on the dimension) such that 

sup P y (S(H(dE k )) ET)< (3.4) 

y&(0,kh(n))\B(Lh(n)) P d 

(see Lemma 5.1 of pQ). Set c n = \M(Lh(n))\, and note that c n < c(Lh(n)) d for some constant 
c. From (13 .3p we have 

W kh{n) (N E , 71 + M fcfe(n) (B(0, fc/i(n))\B(M(n)), T) > M fcMn) (iV I , T) - c n . (3.5) 

We will use Lemma El with L n = M fch(n) (B(0, kh(n))\M(Lh(n)), T) and we note that (HI) 
is ensured by f)3.4p . Let us define 

MT) = E [M kh{n) (N 1 ,71] - £ [M fch(n) (l(0,A:%))\B r ( J L/ l (n)),7l] . (3.6) 

We consider the event that S k is not covered, and use the bound 

1 



P {S k not covered) <P ( 3T C S fc : W fcfc(n) (iV E , T) < -/i(T) 



+ P (s fc not covered, VT C S fc : W feft(ri) (^ I ,T) > ^MT)J ■ 



(3.7) 

In the next Sections, we compute /i(T), and estimate the probabilities of events on the right 
hand side of (13.71) . 

3.1.1 Mean number of explorers crossing a tile 

When T is the tile of a cell C, which has side h(n), and belongs to shell S k C 3(0, n), at a 
distance Ah(n) from B(0, n), then we show that for some positive constants {cd, d > 2} 

MT) > c d A/i(n) d . (3.8) 

The inequality in (I3.8P follows as in (TJ-Section 4.2 and relies on Proposition 12.71 Note 
that d£BD ensures (H2). 

3.1.2 W kh (n)(N I ,T) is unlikely to be small 

Like in (4.17) and (4.18) of Section 4.2 of [I] that there are constants Cd such that 

e ^w^^!* 1 0.9) 

By Lemma l2~4"t since for A large enough we have /x(T) > 3cL d h(n) d > 3c„, there are positive 
constants {c^, d > 2} such that 

P f w f /v 1T T\ <r ^ „m^ <r / exp (-^Ah 2 (n) + X 2 c' 2 h 2 (n) log(n)) for d = 2 , 
P [W kh(n) (N 1 , 7) < -MT) j < | exp (_ A/ ^ (n) + AVXW) for d > 3 . 

(3.10) 
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Thus, after optimizing over A, we get 

P 6k G S fc : ^ (n) (iV 2o, T(*)) < -KT) )<<. d f K ^ hd{n) 



1 \ \ n exp — f-r-, — rx tor d = 2 

../'T^ I ^ J V 4c' 2 log(n) / 



n exp ( — d j , ) for d > 3 . 



(3.11) 

and the event {Wkh(n)(N $-o,T) < §£t(T)} has a probability that decreases, for A large 
enough, faster than any given power of n. 

3.1.3 C is likely to be covered when Wkh(n){N fl , T) is large 

We consider here the event {VT C S, Wkh(n)(N 1 , T) > §^4A d (n)}. Consider shell 5^ at a 
distance A/i(n) from 9B(0, n). Since <Sfc is the union of B(z, h{ri)) when z G Lemma 1X731 
implies, when d — 2, that 



/z 2 (n) 
log(n) 



P (S fe G* A(iV) and W kh(n) (N 1 ,T) > ^(T) for all Tj < \S k \ exp (-k 2 kA 

(3.12) 

For any given power of n we obtain a negligible bound when /i(n) = log(n) and A is large 
enough. When d > 3, then we have 

P (s k G* A(N) and W fcfc(f0 (iV I , T) > ^(T) for all Tj < \S k \ exp (-K d KAh 2 {n)) . 

(3-13) 

For any given power of n we obtain a negligible bound when h 2 (n) = log(n) and A is large 
enough. 



4 Outer error 

In this Section, we prove the outer error estimate ( 11. 5p . This is a consequence of our inner 
error estimates, of Lemma H.5[ combined with coupling with a flashing process of pp. When 
dimension d = 2, and for A large to be chosen later, we decompose the event {5o{n) > 
Alog(n)}, as 

{5 (n) >Alog(n)} = \J{5 (n) G [Alog(n) + z - 1, Alog(n) + z[} . (4.1) 

i>l 

In dimension d > 3, ^/log(ra) replaces log(n) in (14. ip . Note that the index z is at most of order 
n d . Now, we fix i > 1, and we set 3h{n) = A^log^) + z in d > 3, and 3h(n) = A\og(n) + z 
in d = 2. We consider the event {So{n) G [3/i(n) — l,3h(n)[}. We define also, 

£ = 1(0, 7i + 3/i(n))\B(0, n + 3/i(n) - 1). 

Note now that 

P{5 {n) G [3h{n) - l,3h{n)[) < P ((J {z G A(iV), 5 {n) = \\z\\ -n}j . (4.2) 
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For z G E, and in view of Lemma [1.5[ we define 

G(z) = {z G A(N),8 (n) = \\z\\ -n, \A(N) n B(z, h(n))\ > Ph d (n)} . (4.3) 

To prove that P(z G A(N), 8o(n) — II z II ~ n ) is smaller, for A large enough, than any given 
power of n we further split the event into two pieces. 

P(ze A(N), 8 (n) = \\z\\ -n)< P(G(z))+P(z G A(N), \A(N) n M(z, h(n))\ < f3h d (n)) . 

(4.4) 

The second term on the right hand side of (I4.4p is dealt with using Lemma [L~5l We deal now 
withG(z). Note that under {8o(n) G [3h(n) — 1, 3h(n)[}, no explorer escapes B(0, n+3h(n)). 
Thus, on G(z), there are at least (3h d (n) explorers which settle on M(z, h(n)) before exiting 
B(0, n + 3h(n)). We now express the event G(z) in term of flashing explorers, as introduced 
in p. 

4.1 On a flashing process 

We refer the reader to Section 3.1 of pp for a definition of flashing processes. Here, we 
partition 7L d into shells encaging B(0, n), with for k > 

S k = 1(0, n + 2(k + l)/i(n))\B(0, n + 2kh(n)). 

Also, for > 0, let = 9B(0,n + (2A; + l)h(n)). We now consider the flashing process. 
Explorers behave like internal DLA explorers, as long as they stay in B(0,n). After exiting 
B(0, n) they do not flash until their hitting of So, and behave like flashing explorers as defined 
in Section 3.1 of [Tj. In shells {Sk, k > 0}, cells and tiles have the meaning given in Section 
4 of pp. The key features the reader has to keep in mind are as follows. 

• If a flashing explorer is unsettled up to time if (£&), then after time if (£&), it probes 
one site distributed almost uniformly over the cell centered at S(H(Ek)), and settles 
if the site is unoccupied. 

• When an explorer leaves the cell centered on S(H (£&))> h cannot afterward settle in 
Sk, but perform a simple random walk, independent of other explorers, until it hits 
Efc + i. Thus, if we know that an explorer has reached at time t a site of B(0, n + (2k + 
l)h)\M(0,n + 2kh), then it performs after time t a simple random walk, independent 
of its surroundings, until it reaches E*.. 

• We can build the internal DLA cluster, A(N), and the flashing cluster A*(N) using 
the same trajectories S\, . . . ,Sn such that 

TV N 

A(N) = \J{Si(T(i))}, and A*(N) = \J{Si(I*(i))}> (4-5) 

i=l i=l 

and for all i = 1, . . . , N, T*(i) > T(i). This last property, at the heart of our coupling 
argument between the flashing process and the original DLA, is fundamental. It implies 
that if a DLA explorer has crossed a site before settling, then the corresponding flashing 
explorer has also crossed the site before settling. 

Before introducing more notation, let us explain the simple idea behind our estimate. 
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Heuristics Using the representation (14.51) . event G(z) for zGE implies that at least f3h d (n) 
flashing explorers hit M(z, h(n)) before exiting B(0, n+3h(n)). Consider these explorers after 
the moment they enter B(z, h(n)) C Si for the first time. They are behaving as independent 
random walks until they hit Si. Now, a fraction must hit Si on M(z, 2h(n)) fl Si. We show 
that this latter event has a probability we can estimate through the approach of pQ. 

For simplicity, let us call R\ the radius of Si, that is R\ = n + 3h(n). Recall that for 
A C 1(0, Rt) U (9B(0, Ri), we call W Rl (N I , A) the number of flashing explorers which hit 
A before (or as) they hit Si. In this section, the initial configuration is always N lo, and 
we omit this coordinate in W Rl to simplify notation. Under our coupling ( 14.51) . we have 

G(z) C {W Rl (B(z,h(n))) > (3h d (n)} . (4.6) 

Let z' be the closest site of S x to the line (0, z), and note that \\z — z'\\ < 1. Note that a 
fraction of the W Rl (M(z, h(n))) independent random walks in B(z, h(n)) fl B(0, n + 3h(n)), 
must hit Si in a neighborhood of z'. Indeed, first note that since z' G Si, we have 

|3©(*',2/i(n))nB(0,n + 3/i(7i))| > -\dM(z', 2h(n))\. (4.7) 

Now, for any y G dM(z, h(n)), a random walk starting on y, exits M(z' ,2h(n)) on any site 
of <9B(z', 2h(n)) with a probability proportional to (2h(n)) 1 ^ d . Thus, there is a positive 
constant p such that 

inf P y (S (H(dB(z>, 2h{n)))) G <S 2 ) > p. (4.8) 

y£oB(z,h(n)) 

In other words, each flashing explorer stopped on <9B(z, h{n)) before hitting S x has a prob- 
ability at least p to exit Si from Si nB(z, 2h(n)). Thus, there is a positive constant J, such 
that for any large enough integer k 

p(w Rl (M(z',2h(n)))<^k | W Rl (B(z,h(n)))>k) < exp (-Ik). (4.9) 

From (14.6 j) . we have 

\jG(z)c |J {w Rl (W,2h(n)))>^h d (n)} 

U |J {^jB(/,2/ i H))<|/3/ i d (n)andIU i?1 (B(^,/iH))>/3^H}. 

(4.10) 

Let us now define, for any a > 

F(a) = [J {W Rl (M(z,2h(n))) > ah d (n)} . (4.11) 

Thus, from (14. 10p and ( 14. 9p . and for some constant C > 
P( \jG(z)) < P (W R1 (M(z',2h(n))) > P -Ph d (n)) 

+ |Si| sup P (w Rl {M{z',2h{n))) < ^h d {n)\W Rl (M(z,h(n))) > f3h d (n)) (4.12) 
z'eEi V 2 / 

< P + Cn d - X exp (-I^ph d (n)) . 
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It remains to show that for any fixed a, we can find A (defining h(n)) such that P(F(a)) is 
smaller than any given power of n. 

4.2 Estimating P(F(a)) 

Note that by definition of 5/(n), for T C Ei(= dM(0,n + 3h(n)), W Rl {T) satisfies the 
inequality 

W Rl {T) + M Rl {M{0,n- 5r(n)),T) < M Rl {N 1 ,T). (4.13) 
Thus, for some large constant a d to be chosen later, we have 

\ {n) < ad ^ (w Rl (T) + M Rl (®(0,n-a d ^))^) < M Rl (N 1 ,T). (4.14) 
Inequality ( 14. 13j) puts us in the setting of Lemma 12.51 Thus, we first need to compute 



ji(X) = E[M Rl (N 1 ,T)]-E 



M Rl (B(0,n-a d ^),T) 



(4.15) 



Following the same computations as in Section 4.3 of [T], we have for some constants K,c > 
0. 

m < K (a^n^ x < ^h\n). (4.16) 

Secondly, note that as in Section 4.3 of p], we have that for constants {cd, d > 2} 

£ y d il (4,7, 

In optimizing over A in (12. 5p . we find for (other) constants {cd, d > 2} 

P^r: ^m^-)),i'(«-)>^) + -'{-(^j 

(4-18) 

We conclude using the fact, for ad large enough, the first term of the sum in the right- 
hand-side of ( 14.181) is smaller than any given power of n. This was proved in Section [3] for 
the original internal DLA and the same proof can be adapted for the flashing process we 
consider here. The only difference is that we need a stronger version of Lemma 11.31 where 
P(B(0,i?) (jL A(rf)) is replaced by P(B(0, R) (jL A aR (n)) for some large a (this stronger 
version of the Lemma is actually what we prove in the Appendix). Indeed, we can use 
Lemma 11.31 in the context of our flashing process by considering only explorers that do not 
exit 1(0, n). 



APPENDIX 
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A Proof of Lemma 11.3 



We fix n, a configuration of AR d explorers in B(0, R/2), and we choose z G B(0, R). Then, 
P(B(0,R)<£A( v ))<P(M(0,R)<£A qR ( v ))< P(W aR ( V ,z) = 0) (A.l) 

z&(0,R) 

for any a > 1 (in the sequel a will have to be taken large enough). Let L be a large 
positive real to be fixed later, and let ( be the configuration with one explorer on each site 
of 1(0, aR)\M(z, L). We have 

W aR (r),z)+M aR ((,z) > M aR {r],z) - \M(z,L)\. (A.2) 

Note that W aR (r),z) and M aR ((,z) are independent: we are in the setting of Lemma [2T4l 
Assume for a moment that conditions (HI) and (H2) hold, and in addition 

E[M aR (r),z)]-E[M aR ((,z)}>max\3\M(z,L)\, £ P y (H z < H aR ) 2 J. (A.3) 



Then, we have 

P (W aR ( V , z)=0)< exp (-C (E [M QR (^, z)] - E [M aR ((, z)})) . (A.4) 

We next consider separately the case d > 3 and the case d = 2, estimate the expectation of 
AM*?, z) - M aR ((, z), and show (pO}. 



A.l The case d > 3 

We show in this Section that for some Kd > 0, and A large enough, 

E[M aR ( V ,z) -M aR ((,z)} > ^AR 2 » 3\M{z,L)\. (A.5) 

The proof is based on the following classical estimates. There are Oi,a2 positive constants 
such that for any y,z G Z d 

; 7 Ql ,, d _ 2 < Py(H z < oo) < Z 2 _ 2 . (A.6) 
1 + ||y — ^ 1 + 1 1 2/ — z \\ 

Note first that when L is large enough, (HI) holds. Indeed, 

sup P y {H z < H aR ) < sup P y {H z < oo) < - °L_ 2 < - — -, with k > 1. 

y:||z— y||>Z y:||«-y||>i 1 + K 

(A.7) 

We now estimate the mean number of explorers hitting z. 

E[M aR ( v , z)\ - E[M aR ((, z)] = Y, v{y)Py {Hz < H aR ) - ^ P y (H z < H aR ) 

y&(0,R/2) y eM(0,aR)\M(z,L) 

> v(y)Py(H z <H aR )- Yl Py(H z <oo). 

(A.8) 
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Note that for y e 1(0, R/2), we have 

P y (H z < H aR ) =P y {H z < oo) - E y [ 1 H 

a 

> 

~l + \\y- z\\ a -' z 



id— 2 y 



a>2 



l+\\S(H aR )-z 



Id— 2 



> inf 



1 1 • ■■> SU P , , i i • ■> ■ 

yeB(0,i*/2) 1 + ||y - Z\\ d ~ 2 y edM(0,aR) 1 + \\V - A\ 

Now, for a constant a which depends only on a%, a,2, there is k > such that 

inf p (H z < H aR ) > 

y&(0,R/2) 9 R d ~ 2 

Now, using flA.lOp in flA.8j) . we have a constant c such that 



E[M aR { V ,z)-M aR {(,z)\ >AR G 



R d ~ 



a L 



y:\\y—z\\<aR 

>kAR 2 - ca 2 {aR) 2 . 



\y-z 



Id- 2 



When A is chosen large enough, we obtain flA.5[) . 

Finally, there are constants {C^, d > 3} such that for any z G B(0, R) 



C^otR for d = 3 

^ P y (H z < H aR f < I C 4 hg(aR) for d = 4 

3/GB(0,aii) I Cd for C? > 5 



Thus, hypothesis (IA.3[) holds 



(A.9) 



(A.10) 



(All) 



(A.12) 



A. 2 The case d = 2 

We still have 

G aR (y,z) G aR (z,y) 



P y (H z < H aR ) 



G aR (z,z) G aR (z,Zj 



, and G aR (z,y) = E z [a(S(H aR ),y)] - a(z,y), 



where the potential kernel a(., .) replaces Green's function. Note that for < ||,z|| +R 
we have two positive constants K2 and K' 2 such that 

K' 2 \og{2aR) > G B{Zt2aR) > G aR (z,z) > G B(z>R) {z,z) > K 2 \og(R), 

by Proposition 1.6.6 of Lawler [I]. To estimate G aR (z, y), we use Theorem 4.4.4. of |6j 
establishes that for z 7^ 0, (with 7 the Euler constant) 



(A.13) 

< aR, 



|a(0,z)--log(|H 



2 7 + log(8) 



< 



7T 



7T 



Thus, for y e B(0, aR), and < ||z|| < R, and y ^ z 



Gan(z,y) E 

7T 



log 



[ \\SjHq 



R ~ z \ 



\y - A 



<2Ky. 



(A. 14) 
which 

(A.15) 
(A.16) 



Fluctuations for internal DLA 



15 



When y G B(0, R/2) we get 



G aR {z, y)>- log (2(a - l)/3) - 2K g 



71 



(A.17) 



We choose a large enough so that for some constant Ci, we have, for all y in B(0, R/2), 



G a R(z, 

y) > Cl 

Formulas (1A.13I) . (1A.14I) and (1A.18I) together imply that 



(A.18) 



E[M aR (rj, z)\ = Yl v(y)P y (H z < H aR ) > — 

y eM(0,R/2) 



K' 2 log(2aR) 



C 1 AR 2 



y&(0,R/2) 



Using Lemma 3 of j5], we have, for some positive constant C*2, 



K' 2 \og{2aRY 
(A.19) 



E[M aR {(,z)\ < E[M aR (M(0,aR),z)} < 



C 2 {aR) 2 
log(R) ' 



(A.20) 



We need now to choose L to have (HI) satisfied. Note that for y ^ z, ( 1A.16I) and ( 1A.14I) 
yields 



Py(H z < H aR ) < 



-E 



2 l0 J\\S(H aR )-, 
7i \ \\y — z\\ ' 



K 2 \og(R) 

If || 2 — y\\ > Rj log(i?), we obtain, for some constant C 3 

<7 3 log((a + l)logOR)) 



Py(H z < H aR ) < 



log(R) 



(A.21) 



(A.22) 



When R is large enough, we have that (HI) holds for L = R/\og(R). Note that |B(0, L)\ is 
of order R 2 / log(R) 2 and is much smaller than R 2 / log(R). 

Finally we need to control the sum of second moments. Simply note that, from flA.201) . 



Py(H z < H aR ) < E[M aR (C z)] < 



C 2 a 2 R 2 
\og(R) 



a R )\M{0,L) 



(A.23) 



B Proof of Lemma 11.6 



We will choose an h such that R/2h is a positive integer. We divide S = B(0, 2R) \ B(0, R) 
into R/2h concentric shells of height 2h. For k — 1, . . . , R/2h define 



,2R- 2(k - l)h)\M(0,2R- 2kh), and E fe := dB(0, 2R - (2k - l)h). (B.l) 



Also, we set S = B(0, 2R) C . Then, we start on z G <9B(0, 2R) a flashing explorer associated 
with this partition with an explored region V . The flashing setting is much simpler than 
the one introduced in Section 3.1 of pQ. There is an underlying simple random walk, say 
S*, and each shell <Si,c>2, ... is associated with a flashing site. These flashing sites, say 
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{Zk,0 < k < 2R/h} are obtained as follows. We set Zq = z, and for k > 1 we draw a 
continuous random variable R k on [0, h] with density in r G [0, h] i— >■ dr d ~ l /h d : the flashing 
site Zfc is the exit site from B(S , *(if(S/ c )), after time Then, the explorer settles 

on the first flashing site in 5 \ V. The purpose of the flashing construction is that (i) 
the flashing site is distributed almost uniformly inside the ball M(S*(H(Ek)), h), and (ii) 
P Z (H (B(0, R)) < H(V C )) is bounded above by the probability that the explorer crosses S. 

For a small to be chosen later, we say that y G has a dense neighborhood if |B(y, /i) fl 
V| > and we call Dj, their set. There is k > such that knowing that S* has crossed 
Di, . . . , 

• if S*(H(Ek)) Dk, then the probability that S* does not settle in Sj- is smaller than 

K/3, 

• the probability that S*(H(Ek)) G -Dfc is smaller than K\D k \/h 1 (see Lemma 5 of [3]) 
uniformly over the of position the previous flashing site (in <Sfe_i or, exceptionally, on 
the border of Sfc_i). 

Now, the flashing explorer has crossed the annulus S if Z^ G V for all k > 1. In other 
words, 

R/2h 

{#(1(0, R)) < H(V C )} C f) {Z k G V} . (B.2) 

k=l 

By successive conditioning, we obtain 

/ R/2h \ R/2h . 

Pz [ n {Zk ev} )- 1\ v + w) ■ (b,3) 

By the arithmetic-geometric inequality, and ( 1B.2I) . we obtain 

/ R ^\n\\ R/2h 

Note that each y G D k satisfies |B(y, h) D V| > but each site in B(y, /i) fl V is in the 
neighborhood of at most h d ~ l sites of D k . Thus, for some k' 

R ^m^ < . 1 ^ \D k \ „ 2 K '\V\ 

We choose now /3 such that 4/t/3 < 1, and we choose the smallest h such that R/2h is a 
positive integer and 

f /2k/I VI \ ^ 
^>max ho, } . (B.6) 



This adds a constraint on |V|: 



f5 2 R 
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Instead of including (1B.7j) as a condition of our Lemma, we find it more convenient to note 



that the probability we estimate is always less than 1, so that we deal with the case where 
( IB. 71) is violated with the constant of 1 11.111) . 



C Proof of Lemma 11.5 



Recall that and are the constants appearing in Lemma 11.61 We define a positive 
constant 

7=max ( 1; (^) d y (ai) 

Choose now P > such that 4 d (3>y < 1, and h = R/A > 1. Note that 

tM < iPR d < K- (c.2) 

We build now, by induction, a random subdivision of B(z, i?) into shells of heights h , hi, ... , 
in which respectively N ,Ni,... explorers of A(rj) have settled. We emphasize that the 
randomness comes from A(rj), and that the event {0 G ^4(^)} imposes to have iVj > |_^«J , for 
z > 0. Assume that hi, . . . ,hk have been defined such that 

/a* > 1 and V/i;<-. ( C - 3 ) 

i=l 

We define /if+i = ^N k < j\r)\, and, by (10.2|) we have /i^+i < h . Note also that h k+ i > 1. 
Indeed, necessarily iV fc > [^-fcj, so that > 7[^fcJ > L^fcJ- Since mm(hi, . . . , h k+1 ) > 1, 
the number of steps before we violate (1C.3j) . say L, is finite. Obviously L < R. Note that 



since hi < h 

a <J2 h ^ h L+Y, h ^j + ^ ( c - 4 ) 



8=1 1=1 

Thus, we define 



h L+ i = R- [J2hi >0. (C.5) 



,i=0 



For any choice of integers l,n , . . . ,ri[, the event {X = Z, iVo = no, ■ ■ ■ , N L — rii} implies that 

ni H \-rii explorers have crossed a shell B(z, R)\B(z, R — h ) by stepping on at most n 

explorers settled in it, that n 2 + ■ — h n L explorers have crossed shell B(z, R — h )\B(z, R — 
ho — hi) with rti explorers settled in it, and so on and so forth. Using Lemma [1.61 the fact 
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that rii < f3R d , I < R, and the notation 5 = ^-j-, we reach the following estimate. 

P(0 G A(rj)) 

< P(L = l,N = n ,...,N L = m) 



l<R,no,ni,...,ni<\rj\ 
\/i,ni>\hi\ 



< R((3R 



d\R+l 



sup 

l<R,riQ,ni,...,ni<\r] 
Vi,rii>{hi\ 



id 

"o 
n 



+ ...+ 



(C.6) 



Now, note that by the arithmetic-geometric inequality, for 1 < % < I (and using hi < ho) 



M, 

n 



h d 

n i-l 

Ui-l 



> I — X • • • X 

n 



n i-i 



h d 



ni-i 



-7 



Ui-l 
S/i 



(C.7) 



Thus, from flC.6j) and flC.7j) . we have 

P(0 G A(rj)) < R((3R d ) R+1 



max exp | — a d 

l<R,no,ni,...,ni<\rj\ > 
V«,n 4 >L/!.jJ 



i=l 



< P(/3P 



max 

l<R,no,ni,...,ni<\r]\ 
\/i,ni>[hi} 



exp -— J2 ih i+i 



(C.8) 



7 



Since /ii < P/4, note that we have h 2 + ■ ■ ■ + h L > R/A by (1C4[) . By Holder's inequality, 
note that for constants {q, d > 2} 



L-l 



i+l 



> 



(Ef=i 1 



8=1 



1 

2-n=l jl/(d-l 



d-1 



> 



C 2 



^ - /, '" > lor r/ 2 



>c 2 



(L) - ^log(fl) 



c d7 ^2 > c d R 2 for d > 3 . 



(C.9) 



This concludes the proof. 
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